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Rotating magnetohydrodynamics 


By X. SHAN anD D. MONTGOMERY 


Department of Physics and Astronomy, Dartmouth College, Hanover, 
New Hampshire 03755-3528, U.S.A. 


(Received 25 April 1994) 


It is shown that rotation of a periodic column of magnetohydrodynamic fluid 
can be stabilized against current-driven instabilities by rotation. The rotation 
is considered as resulting from a slight departure from overall charge neutrality, 
and is braked by wall friction (loss of angular momentum to the cylinder wall). 
Both effects are modelled phenomenologically in the equation of motion, and 
the emphasis is on their bulk dynamic effects rather their microscopic kinetic 
theory origins. Much as rotation is known to stabilize a Navier-Stokes fluid 
against transverse displacements, it is shown that sufficient rotation will 
suppress the helical vortices and helical deformations of the current channel 
that are known to result when the axial current in the column exceeds its 
stability threshold. 


1. Introduction 


When a Navier-Stokes fluid is ‘spun up’, or rotated at a fixed angular 
velocity, the fluid elements acquire an effective elastic binding to their steady- 
state positions in the two dimensions perpendicular to the axis of rotation. 
Small oscillations thereby become possible that were not possible for a non- 
rotating fluid (see e.g. Batchelor 1967; Greenspan 1968; Acheson 1990). There 
is stabilization against small displacements of the fluid elements by virtue of the 
rotation. The purpose of this paper is to present evidence that a related 
stabilization through rotation may occur in magnetohydrodynamics (MHD). 

Usually, rotation of a Navier-Stokes fluid occurs because of the rotation of 
a material surface with which it is in contact. The surface might be a cylindrical 
container, for example, or perhaps the bottom of the ocean. Magnetofluids may 
acquire rotation velocities through less direct means. A small departure from 
charge neutrality, for instance, may result in a radially outward or inward 
electric field from Gauss’s law, in cylindrical geometry. If an externally 
supported axial magnetic field is present, the resulting E x B particle drift can 
produce an azimuthal fluid velocity. Rotation of plasmas in ‘Q machines’ has 
been so interpreted since the 1960s. Chen & Montgomery (1993) have recently 
treated the case of a slightly non-neutral magnetized plasma in an axial 
magnetic field in thermal equilibrium, and have identified a rigidly rotating 
axisymmetric state as the thermal equilibrium one in the presence of perfectly 
smooth impenetrable walls in a periodic cylinder geometry. The charge excess 
or deficiency was shown to be concentrated in a boundary layer near the 
cylinder wall. However, the necessarily complex microscopic interaction of real 
plasmas with their boundaries (limiters or divertors, for example, in the case of 
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tokamak confinement) will make the rigidly rotating, or ‘spun-up’, profile an 
inadequate description. 

Both the plasma/wall interaction and the radial distribution of a small net 
charge in a magnetized plasma are inherently kinetic processes, and their 
inclusion in a self-contained MHD description will necessarily have a somewhat 
ad hoc phenomenological quality. Nevertheless, it seems likely, because of their 
complexity, that such a phenomenological inclusion in MHD is the only 
possibility for assessing their effects on fully general, time-dependent, nonlinear 
behaviour. 

An experimental situation that led to the investigations reported here (but 
otherwise not connected to it in detail) is the so-called ‘L-to-H transition’ in 
tokamaks. ‘L’ stands for ‘low’ and ‘H’ for ‘high’, and these refer to the degree 
of confinement and quiescence in the plasma. An apparently spontaneous 
transition was sometimes observed in tokamaks (see e.g. ASDEX team 1989) 
to a regime (the ‘H mode’) with superior confinement properties. The 
desirability of ‘H-mode’ operation has generated a large literature of further 
experiments and proposed explanations (an on-line literature search of the 
references since 1988 returns about 150 publications; a useful review by Burrell 
et al. (1992) contains over 80 citations). The phenomenon is complex, and 
exhibits a morphology that varies from machine to machine. A single detailed 
explanation encompassing all of the observed phenomena seems increasingly 
unlikely. We focus here not on the experimental L/H transition as such, but 
rather on an MHD problem that one subset of the L/H measurements has 
suggested : stabilization through rotation. 

An early observation of the L/H transition was that it was often accompanied 
by the appearance of a flow velocity, often but not always in the azimuthal 
(‘poloidal’) direction. A ‘shear’ was often inferred in the velocity field (i.e. a 
departure from a rigid rotation profile) from radial electric field measurements 
and their interpretation in terms of an E x B drift. Taylor et al. (1990; see also 
Tynan et al. 1992) showed explicitly how an electrostatically produced poloidal 
rotation velocity could trigger the H-mode transition. Theoretical inter- 
pretations of this observation have been diverse and full of puzzles and 
controversy. 

The focus here is on what we believe to be a generic MHD result suggested 
by some of these observations: that stabilization of a current-carrying 
magnetofluid can result from poloidal rotation (of whatever origin), in much the 
same way that rotation is known to stabilize a Navier-Stokes fluid. The method 
is numerical solution of the MHD equations, suitably modified to include 
rotation and wall friction. A fully three-dimensional MHD spectral code 
(previously described in detail in Shan, Montgomery & Chen 1991; Shan & 
Montgomery 1993a, b; Montgomery 1993) can be readily adapted to include the 
rotation and wall friction. The two modifications to MHD are described in §2, 
which also summarizes previous results obtained from the code for the non- 
rotating case. Section 3 is a brief summary of the numerical method, and can 
be skipped by a reader who is interested only in the results. Essentially no new 
numerical complications result from the modifications necessary in the 
numerics. Section 4 is devoted largely to the changes, from the previously 
published results, to which the rotation and wall friction lead. Section 5 is a 
brief summary, and identifies possible future directions. 
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2. Rotating MHD with wall friction 


We write the MHD equation of motion in a standard set of dimensionless 
(‘Alfvénic’) units. Only the first and last terms on the right-hand side are in any 
way non-standard: 


Ca ee 


VoL; = Quy — 
z : +jx B—Vp+vV*v—A(r)v. (1) 


Here v = v(x, t) is the fluid velocity, assumed solenoidal; the mass density of the 
magnetofluid is assumed to remain spatially uniform. The mechanical pressure 
is p, the magnetic field is B, j = V xB is the current density, and v is the 
kinematic viscosity. The function v, is a poloidal (azimuthal) rotation velocity, 
and is assumed to be given as a function of position. 7 is a spin-up time, also 
assumed to be given. For vo, we have typically chosen v, = fro, €, in cylindrical 
co-ordinates (r,¢,z), with w, a uniform axial vorticity. 7 is also assumed to be 
a uniform constant. We have experimented with assigning radial variation to 
v, and 7, but without a noticeably significant change in the results. A = A(r) is 
a phenomenological friction coefficient, which only becomes significant in 
magnitude near the radius of the cylinder, which lies at r = a. This last term in 
(1) models what are complicated, situation-dependent kinetic processes 
involving loss of momentum and angular momentum to whatever mechanical 
element (limiter, divertor) is being idealized at r = a. Without the — A(r) v term, 
the first term on the right-hand side of (1) would by itself have the effect of 
‘spinning up’ the magnetofluid to the azimuthal velocity v, in a time of order 
T. It is understood to model the effects of non-neutrality in a way that 
consciously oversimplifies the kinetic processes involved in the interests of 
having a tractable representation of their effects in bulk MHD. The wall friction 
term frustrates the attempt to achieve rigid rotation, and imposes a ‘shear’, or 
radial gradient in the axial vorticity, in the region where the first and last terms 
on the right-hand side of (1) compete. 
The magnetic field is advanced by Faraday’s law, 
ôB 


gT VE (2) 


where E is the electric field. Ohm’s law, in the form 
E+vxB=nj, (3) 


is assumed to apply. The magnetic diffusivity 7 is assumed to be a given 
function of the radius r that increases steeply towards r = a, reflecting the 
strong dependence of electrical conductivity on temperature. The viscosity v, 
however, is treated as a uniform constant. 

The geometry considered is a straight cylinder of length L,, with periodically 
identified ends — the lone concession to the toroidal geometry in which the L/H 
transition has typically been studied. An axial magnetic flux with mean 
magnetic field B, is assumed to be supported externally. An applied electric 
field E = Ẹ,ê, is assumed to drive a current in the axial direction. In the 
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absence of any mechanical activity (v = 0) or spin-up (v, = 0), the steady-state 
current density is j = Z,€,/7 and the azimuthal (poloidal) magnetic field is 


Blr) = (27r) [2 rdrd¢, 


where the surface integral is over the area of a circular cross-section of radius 
a. The pinch ratio is defined as 0, = By(a)/By, and the Hartmann number 
(evaluated at r = 0) is H, = aC ,/(nv)?, where C4 is the toroidal Alfvén speed (H, 
is just B,/(yv)! in these dimensionless units). It has been shown previously 
(Montgomery 1992; Shan & Montgomery 1993a) that ©, and H, govern the 
linear stability threshold of any axisymmetric, non-rotating state. ©, may also 
be expressed in terms of the ‘safety factor’ g(a) evaluated at r= a: 


27a 


294(4) 

The wall at r =a is imagined as a rigid perfect conductor, mechanically 
impenetrable, and coated with a thin layer of insulating dielectric at r = a. The 
boundary conditions implied by this are that at r=a, j.ñ =0, B.i=0, 
v. = 0, where ñ = 6, is the unit normal to the wall. A second mechanical 
boundary condition may be chosen, and a classical fluid-dynamical perspective 
might suggest either v x ñ = 0 (‘no slip’), or ñ x ø. Â = 0, where ø is the viscous 
stress tensor (a ‘stress-free’ wall). With the numerical method used, both these 
boundary conditions turn out to be quite difficult to implement, and probably 
represent no better the complex plasma/wall interaction than the one we do 
choose: w. Â = 0, where w = V x v is the vorticity. The condition w.f = 0 is 
implied by, but does not imply, no-slip boundary conditions. The wall drag will 
be simply accounted for by the term — A(r) v in (1). A(r) becomes significant only 
near r = q. 

Previous publications (Montgomery, Phillips & Theobald 1989; Chen, Shan 
& Montgomery 1990; Agim & Montgomery 1991 ; Shan et al. 1991 ; Montgomery 
1992, 1993; Shan & Montgomery 1993a,b) have explored in some detail the 
dynamical events that occur as Æ, and thus ©, are raised from zero. A linear 
stability boundary in the (H,,@,) plane can be sharply identified for a given 
y(r). The linear stability analysis can be given analytically for a uniform 7, but 
requires numerical computation for a spatially variable one (Shan & 
Montgomery 1993). As long as ©, remains below the value where any normal 
mode grows temporally, arbitrary initial perturbations damp away owing to 
resistivity and viscosity, and the magnetofluid returns to its axisymmetric, 
zero-flow, steady state. 

When Q, is raised just above its (calculable) stability threshold, the steady 
state bifurcates. The preferred branch (as revealed by dynamical computation) 
is no longer axisymmetric, but helically symmetric. The MHD fields are 
additive combinations of the axisymmetric steady-state fields and their linearly 
unstable eigenmodes, plus small higher-order contributions from their inter- 
actions. In this regime an arbitrarily perturbed initial state will relax to a 
time-independent steady state with helical symmetry and a lower energy 
dissipation rate (at fixed Æ) than the corresponding axisymmetric resistive 
equilibrium, which is now unstable. 

Further increases in Æ, lead to other interesting transitions, not yet fully 
mapped out. For example, the poloidal and toroidal mode numbers (m and n) 


o= 
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that characterize the helical state may change (the unstable eigenmodes vary 
as exp (tm @—1i2mnz/L,)). ‘Multi-mode’ states, with more than one m/n ratio 
present, may be excited, and may acquire a persistent quasi-periodic time 
dependence. Eventually, as Æ, is increased further, fully developed MHD 
turbulence sets in, and has been both difficult to characterize and difficult to 
compute. Our numerical methods are not adequate to explore the fully 
turbulent state. 

In summary, as Æ, is raised, an increasing level of mechanical and magnetic 
MHD activity sets in, and at each transition the fields seem less likely to be ones 
that would confine a plasma well or minimize its transport. The analogy with 
hydrodynamic pipe flow, as the pressure drop and Reynolds number are raised, 
is non-trivial, and the observed MHD behaviour is of a piece with hydrodynamic 
scenarios as viscous flows are subjected to steeper pressure gradients. 

The question of interest here is taken to be ‘What possibilities exist for the . 
two new terms in (1) to suppress this MHD activity, and what levels of rotation 
w, are required?’ How does one characterize the rotating state? Easy and 
relatively superficial answers, for the case of uniform 7 and A = 0, can be given 
analytically (Montgomery, Chen & Shan 1991), but for the case of radially 
increasing 7 and non-zero A numerical computation is required. We now give a 
brief summary of the numerical method before turning to the results to which 
it leads. 


3. Numerical method 


The code has been described elsewhere in considerable detail (Shan et al. 
1991; Shan & Montgomery 1993a); this section is a brief description of its 
operation. All fields are expanded in Chandrasekhar—Kendall functions, or 
orthonormal eigenfunctions of the curl. These start from any solution of the 
scalar wave equation 


(V2+A2) yy = 0, (4) 
and are defined, before normalization, by the expression 
A=AVxéy+Vx (V x ey), (5) 


where ê is a constant unit vector. It is easy to show, using (4), that Vx A = AA. 
For the case of the periodic PTa the appropriate choice of ê is ê,, and 


y= Vang = Jm(Ynma”) exp( (img —ikn z). 
Here m = 0, +1, +2, ...; n= 0, +1, +2, ...; k, = 2an/L,; q = 1, 2, ...; and 
Vieng t Bi, = À? mg defines Apmg (both signs allowed). Thus the A,,,,, are indexed 
by three integers, and may be normalized so that 


(Anmg Anima) =$ [Abn As mgt U Ob dz = Sin Sm Sag (6) 


where the volume integral runs over the interior of the periodic cylinder, and 
V = maL, The Ynmq are determined by imposing the boundary condition 
Anma:€r = 0 at r =a, for all m, n, q such that m?+n? > 0. For m =n = 0 the 
Yoog are determined by demanding that f A,,,.é¢,rdrd¢ = 0 when integrated 
over a cross-section of the cylinder. This makes all the Ago, ‘fluxless’, and, since 
one is interested in expanding fields such as j and B that have non-zero axial 
fluxes, the Apoy must be supplemented by a function A,(r) with a finite flux, 
orthogonalized to all the retained Ase by the Schmidt orthogonalization 
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procedure (for details see Shan et al. 1991). All the fields are then expanded in 
the A,,,, and A,(r); for example, 


nma 
B = jy Ag +B,é,+ x fani oe (7) 
nmga 
V =w Atv + Z Er na Anma (8) 
nmg 


The £2.19: Snmq Jo and w are time-dependent complex expansion coefficients, 
and serve as the dynamical variables of the problem. 

The expansions (7) and (8) and their curls are substituted into (1)-(3). Then 
inner products are taken one at a time with the expansion functions, and (6) 
permits ordinary first-order differential equations to be extracted in the usual 
fashion of a Galerkin approximation. These have the generic form 

aX, _ XAR 

ae D Cap X gt BC apy Xp Xyt = (9) 
where X, is an element of a column vector whose elements are all the real and 
imaginary parts of the £3..4, Ehmo Jo and wo. The coupling coefficients C,, and 
Cs, involve integrals of Bessel functions, and are stored as a table. The Æ, stand 
for any external forcing terms that may be present (in this case the effect of 
applied axial electric field). The number of independent X, retained may 
typically be about 2000, and they are advanced according to (9) by a time- 
stepping routine appropriate to ‘stiff’ ODEs. The method converges rapidly 
(i.e. the answers become independent of the further addition of expansion 
functions) for situations not too far above the stability threshold in ©,, but does 
not converge well for fully turbulent cases. 

The time advancement takes place entirely in spectral space. Only when it is 
necessary to construct configuration-space graphics it is required to reconstruct 
the fields as functions of x from (7) and (8). The reason the code converges 
rapidly, when it does, is that the expansion functions are close to the linear 
normal modes of the system. They are exactly the linear normal modes for the 
case of uniform 7, and the m and n indices are still ‘good quantum numbers’ 
even when 7 depends upon r. Often, a numerical code does not have the luxury 
of a set of expansion functions that individually are already quite close to the 
linear eigenfunctions of the physical problem, including the boundary 
conditions. For example, in the case of plane shear flows, the Orr-Sommerfeld 
eigenfunctions are far from the Fourier-Chebyshev expansion functions 
sometimes used to decompose them (see e.g. Canuto et al. 1988). 

The two new terms in (1) are readily accommodated in the numerical scheme, 
and essentially no innovations are required. We turn now to the differences in 
the dynamical behaviour that the new terms induce in the near-threshold 
dynamical behaviour, relative to that behaviour in the absence of rotation and 
wall friction, as previously described. 


4. Rotation damping of helical MHD activity 


In order to compare the rotating MHD results with the non-rotating ones, we 
pick again the resistivity profile used by Carreras et al. (1980) and Hicks et al. 
(1981) in studying ‘tearing-mode’ activity. This was used earlier by us (Shan & 
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B, v No A, E, O An Ta 
Run I 4 107 10-74 40000 139x10 0052 10 09 
Run II 2 2x10 1074 14142 116x10% ©i1 4 07 
Run HI 2 2x10 1074 14142 146x10% O11 10 0.9 


TABLE 1. Parameters of the runs (7, and H, are evaluated on axis) 


Montgomery 1993a, b) in the investigation of the dependence of MHD activity 
upon the on-axis Hartmann number. For the magnetic diffusivity profile, we 


take 
r 2A]I+1/A 
nr) = nfi+(2) | , 
To 


with a = 1-0, 7, = 0°6, aspect ratio L,/2ma = 2:5 and A = 4-0. The steady-state 
current profile is 
jr) = E,/y(r). 


The stability boundary in the plane has previously been determined and 
exhibited (Shan & Montgomery 19936). The safety factor profile 


g(r) = q(0) [i + oT 


To 


has the value q(0) = 47aB,7,/E,L, on axis, is relatively flat out to the 
neighbourhood of r = r, and then rises more abruptly toward the wall r = a. 
For the wall friction coefficient A(r), we must also make a definite choice, and 


we pick the function 
r 2A, )I+1/A, 
A(r) =77 |: +(2) | =r}, 
Tx 


with A, and r, as constants to be specified. This function is almost zero in the 
interior, and becomes significant only at the edge of the cylinder, r > r. The 
viscosity is treated as a uniform constant throughout. (We do not inquire here 
into how the net plasma charge arises or distributes itself.) The reciprocal of the 
spin-up time 77? is allowed to go from zero to a finite value of 0-01 at the time 
it is desired to begin studying the effects of rotation. The angular rotation rate 
lw, may be changed at times of choice when it is desired to study the effect of 
increasing the rotation rate. All times are measured in Alfvén transit times of 
unit distance and unit magnetic field strength. 

To illustrate our results, we have made three runs with the code, hereafter 
called runs I, II and III. Table 1 summarizes the various parameters used in 
these three runs. All these calculations are started with the trivial equilibrium 
v = 0 and j = E, ê,/n(r), plus small-amplitude (of order 10 in amplitude) noise 
randomly distributed amongst all the expansion functions A,,,7. This trivial 
equilibrium is unstable with all three sets of parameters given in Table 1. The 
two new terms in (1) are turned off initially until the helical unstable modes 
grow and saturate. We then study the effects of rotation on the helical MHD 
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FIGURE 1. Contours are intersections of magnetic surfaces (a, c, e) and vorticity field line 
surfaces (b, d, f) with the plane z = 0 at three different times for run I: (a) and (b) are at 
t = 600, after the helical MHD activity has developed and at the instant the rotation is 
turned on; (c) and (d) are at t = 700; (e) and (f) are at t = 1200, at which time the helical 
MHD activity has been almost totally suppressed. Times are in Alfvén transit times of unit 
magnetic field and unit distance (a = 1). 
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activity by turning on the wall friction effect and raising the rotational driving 
term incrementally. 

For run I the parameters are barely on the unstable side of the stability 
boundary of the m/n = 3/1 mode in the (@,,H,) plane. Small initial random 
noise is allowed to amplify, and, as before, the magnetofluid goes into a helical 
steady state with stationary paired helical vortices with m/n = 3/1. Since this 
steady state is purely helically symmetric, well-defined flux surfaces exist for 
the quantities B, j, v and w. Figure 1 (a) is a computer-drawn cross-section of the 
magnetic surfaces that characterize this steady state, and figure 1 (b) is a cross- 
section of the flux surfaces of the steady-state vorticity, both at time t = 600. 
At time ¢ = 600, 77! is allowed to take on the value 0-01, with w = 0-2. By 
t = 700 (figures 1c, d), the departure from axisymmetry in these two quantities 
has been noticeably suppressed, and by ¢ = 1200 (figures 1e, f), a nearly 
axisymmetric state has been achieved. This suppression process of the helical 
deformation can be seen more clearly in the modal energy plots shown in figure 
2, where each line in (a) and (b) is the time history of the kinetic or magnetic 
energy of the helical modes with a specific m/n ratio, i.e. the same helical 
symmetry. From the two contour plots of w, shown in figure 3, at t = 550 and 
1000, it can be seen that a ‘shear layer’, or layer of steep radial gradient in the 
axial vorticity has developed just inside the wall r = a = 1 at t = 1000, at which 
time the axisymmetry has been recovered almost totally. The behaviour 
illustrated in figures 1-3 is typical of what happens if sufficient rotation is 
applied to stabilize the magnetofluid completely against the appearance of the 
helical states. The azimuthal velocity profile at t = 1200 is shown in figure 4. In 
terms of the toroidal Alfvén speed, the amount of rotation applied in this case 
to make the axisymmetry disappear corresponds to a maximum rotation speed 
of 0:03 times the Alfvén speed. Note that without the wall friction term in (1), 
the radial dependence of the azimuthal velocity would be linear. The boundary 
layer begins just above rą = 0-7 here, the choice of which is optional. 

In run I we have demonstrated the suppression of a purely helically 
symmetric stationary mode by the rotation of the magnetofluid. More generally, 
the rotation will make the magnetofluid inside the cylinder more stable against 
any perturbations, even when the helical MHD activity was not stationary. We 
have previously shown that when either ©, or H, is raised to sufficiently higher 
values, the magnetofluid will bifurcate into some more and more complex 
states. We shall demonstrate next in run II that the rotation of the magnetofluid 
has an effect opposite to that of increasing ©, and H,. The parameters for run 
II put it well inside the unstable region of both 2/1 and 3/2 modes in the (0), Hy) 
plane. (This corresponds to what was called run 3 in Shan & Montgomery 
(1993b).) As before, in the initial stage, when both the rotational driving force 
and the wall friction are turned off, the unstable modes grow out of the random 
perturbation and form a multi-mode state with m/n = 2/1 and 3/2 being the 
largest modes. We then turn on the wall friction and increase w, gradually as 
a sequence of step functions. The gradual suppression of the helical 
perturbations is most apparent in the time history plot of the modal kinetic and 
magnetic energies of the helical perturbations shown in figure 5. As the rotation 
rate is increased, the magnetofluid becomes more and more quiescent: from the 
three-dimensional state to two-dimensional helically symmetric state, and 
eventually to the one-dimensional axisymmetric state. Since in states involving 
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Ficure 2. Time histories of the modal kinetic energies (2) and modal magnetic energies 
(b) for run I. The traces are for all modes with given m/n, summed over q. 


many different m/n ratios the existence of any flux surfaces is impossible to 
assert, we show in figure 6 a series of contours of the axial vorticity field w, at 
typical times (a) 300, (6) 900, (c) 1500, (€) 2100, (e) 2700 and (f) 3300. The 
transition to a helically symmetric state (d, e) and finally to an axisymmetric 
equilibrium (f) is again apparent in this series of plots. 
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Figure 3. Contours of constant axial vorticity component in the plane z = 0 at times 
t = 550 (a) and t = 1000 (b) for run I. 
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Fiaure 4. Azimuthal velocity field averaged over a cylindrical surface of radius r, plotted as 
a function of r, for run I at ¢ = 1200. The maximum rotation velocity is about 0-03 of a 
toroidal Alfvén speed. 


The third run would otherwise be identical to run II, except that here we 
employed a thinner wall friction layer by choosing different A, and r,. An 
interesting feature of this run is the sawtooth-like oscillation in the less- 
dominant m/n ratios when the rotation induces a transition from the ‘mixed- 
mode’ state to the helically symmetric state, as can be seen in the modal energy 
plot shown in figure 7. This oscillation, which occurs between about t = 1000 
and ¢ = 3500, is a limit-cycle solution of the dynamical system (9), and has been 
reported in previous calculations (Shan & Montgomery 19936). Still further 
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Figure 5. Time histories for run II of the kinetic modal energies (a) and magnetic modal 
energies (b), with all modes corresponding to a singe m/n combined. For this run, the rotation 
rate is stepped up incrementally from zero, at the instants indicated by the arrows at the top 
of the graph. The values of axial vorticity turned on are printed just after the arrows. 
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Figure 6. Contours of constant axial vorticity in the plane z = 0 at times t 
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Ficure 7. Kinetic modal energies (a) and magnetic modal energies (b) for run III. This run 
is for the same plasma conditions as run II, except that a thinner wall friction layer has been 
induced by choosing a value of rą =0°9 and A, = 10. The rotation is again increase 
incremen tally, with the times indicated by arrows at the top of the graph, and values of the 
axial rotation vorticity w, printed immediately to the right of the arrows. 
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increases in the rotation rate would lead to the total disappearance of all the 
helical contributions. The single-m/n state might be called the ‘MHD L mode’. 
The MHD system oscillates around the helically symmetric equilibrium, in a 
way quite similar to what has been observed in experiments by Zohm et al. 
(1994) during an L-H transition. 


5. Discussion and possible extensions 


We consider it noteworthy that rotation seems to have a generic effect of 
stabilizing current-driven MHD activity near its stability threshold. It is 
perhaps significant that the layer of strong velocity shear, visible for example 
in figure 3(b), has appeared as a consequence of the combined effects of the non- 
charge-neutral rotation and the boundary friction, rather than as a cause of the 
rotation, as has sometimes been suggested. The modelling done for stepping 
outside the MHD framework to get the E x B rotation due to net charge and the 
wall interaction is crude, and it may be a very long time before microscopic 
kinetic theory will advance to the point at which the adequacy of the model 
terms in (1) can be tested. Nevertheless, the experimental message is clear: 
rotation is a good thing from the point of view of bulk MHD stabilization. Since 
experimental techniques are available for inducing rotation (by, for example, 
the introduction of electrons into the column, as in Taylor et al. 1989), it would 
appear to be desirable to study empirically the effects of making confined 
plasmas rotate in a systematic fashion. 

The limitations to uniform density and incompressible flow, the absence of 
Braginskii viscous stress tensors (particularly the ‘ion-parallel’ part) and the 
lack of a local temperature dependence for the conductivity are all severe 
limitations, the lifting of which would be desirable. However, none of them is 
trivial. 
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